Total light absorption in graphene 
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We demonstrate that 100% light absorption can take place in a single patterned sheet of doped 
graphene. General analysis shows that a planar array of small lossy particles exhibits full absorption 
under critical-coupling conditions provided the cross section of each individual particle is comparable 
to the area of the lattice unit-cell. Specifically, arrays of doped graphene nanodisks display full 
absorption when supported on a substrate under total internal reflection, and also when lying on 
a dielectric layer coating a metal. Our results are relevant for infrared light detectors and sources, 
which can be made tunable via electrostatic doping of graphene. 
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Beyond their fundamental interest in quantum theory 
and cosmology, blackbodies have practical relevance in 
light absorption and emission devices (e.g., detectors, 
photovoltaics, and optical sources). Inspired by nature, 
the following two routes have been investigated to achieve 
perfect absorption: (i) diffusion in disordered lossy sur- 
faces, such as black silver, has led to engineered materi- 
als made of carbon nanotubes that exhibit extraordinary 
broadband light absorption [1] ; (ii) alternatively, ordered 
periodic structures, which explain the moth eye effect 
observed in nocturnal insects |2] , have been pioneered by 
experimental and theoretical studies showing total light 
absorption (TLA) in the visible using metallic gratings 
[2H5] . Similar phenomena have been reported at infrared 
(IR) [6] and microwave [7] |8j frequencies. Recently, om- 
nidirectional TLA [5] has been realized using periodic 
surfaces supporting localized plasmon excitations. 

The availability of high-quality graphene as a stable 
material with extraordinary (opto-)electronic properties 
[TUl - fT2"] makes a compelling case for exploring its ability 
to harvest light for potential application to optoelectron- 
ics, with the advantage of being optically tunable via 
electrostatic doping [13]. However, a single sheet of ho- 
mogeneous graphene is poorly absorbing [T3] (about 2.3% 
absorption), so the challenge is to transform it into a per- 
fect absorber, for which we can rely on its power to host 
extremely confined plasmons [HI [16] . 

In this Letter, we show that a single sheet of doped 
graphene, patterned into a periodic array of nanodisks, 
exhibits 100% light absorption. We first discuss the ex- 
tinction cross-section of graphene disks, which can exceed 
by over an order of magnitude their geometrical area |16j . 
These disks therefore belong to the class of absorbing par- 
ticles that can be arranged in periodic planar arrays such 
that their cross section exceeds than the area of the unit 
cell. We further assume here that all diffracted beams are 
evanescent. Under these conditions, we derive a universal 
maximum absorption for any thin layer separating two 
different media. In particular, TLA is predicted when 
the transmission channel is suppressed, for instance un- 



der total internal reflection (TIR) by illumination from 
a prism substrate, and also when the particles are sup- 
ported on a dielectric film on top of a thick metal. We 
show that a patterned graphene sheet displays these ef- 
fects. Full absorption by an atomically thin carbon film 
is thus possible using currently available technology. 

Optical response of graphene nanodisks.- Doped 
graphene has been predicted to support strongly con- 
fined, long-lived plasmons |15j . which produce sharp IR 
resonances in nanodisks tens of times smaller than the 
light wavelength [TB]. Figure Qa) shows an example sim- 
ulated with the boundary element method (BEM) [17] . 
We then fit the polarizability of these small disks to [13] 

a(uj) = — — 5-^, (1) 

where lo p is the plasmon frequency, K is its decay rate, 
and K r is the radiative contribution to k. The extinction 
cross-section a cxt = (Attuj / c)lm{a} obtained from a [Fig. 
[lja), broken curve] is in good agreement with rigorous 
numerical solutions of Maxwell's equations (solid curve). 
We represent the parameters w p , k, and n r in Fig.JIJb-c) 
for different doping levels (quantified by the Fermi en- 
ergy Ep) as a function of disk size. Here, we compute 
the conductivity of graphene in the local-RPA limit [15] 
with an intrinsic relaxation time r = fiEp/evp, where 
vp ~ c/300 is the Fermi velocity and fi — 10, 000 cm 2 /Vs 
is the measured DC mobility [10] (e.g., t s» 10~ 13 s 
for Ep — 0.1 eV). Interestingly, the maximum extinc- 
tion cross-section a^fax ~ (3A 2 /27r) k t J k at the resonant 
wavelength A can exceed by over an order of magnitude 
the geometrical area, although the disks are far from per- 
fect two-level systems (such as an atom without inelastic 
decay channels), for which c™* x ~ 0.48A 2 . This is a con- 
sequence of the fact that k » k,. 

Universal limit to absorption by a thin layer.- In a 
symmetric environment, the reflection coefficient r of a 
thin material layer (much thinner than A) determines the 
transmission coefficient 1 ± r (upper and lower signs are 
for s- and p-polarized light, respectively, as required to 
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FIG. 1: Optical response of doped graphene nanodisks. (a) Characteristic extinction cross-section of a nanodisk (60 nm in 
diameter, and Fermi energy Ef = 0.4 eV), dominated by a pronounced dipolar plasmon and strong near-field enhancement. 
Solid curve: full electrodynamic calculation. Dashed curve: fit to a Lorentzian polarizability. Inset: near-field intensity relative 
to incident intensity for plane-wave illumination with polarization parallel to the disk. (b,c) Size-dependence of the maximum 
extinction cross-section and the polarizability parameters [see Eq. Jy]. 



enforce the symmetry of the scattered field with respect 
to the layer plane; we use this convention from now on). 



Therefore, the maximum absorption (1 



|l±rf) 



is limited to 50% (corresponding to r = =pl/2). This 
value can be increased by considering asymmetric envi- 
ronments, such as the patterned graphene sheet depicted 
in Fig. [fa). 

A simple analysis leads to the maximum possible ab- 
sorption produced by a thin, periodically structured ma- 
terial layer that is infinitesimally close to the interface 
between two media of different permittivities, t\ and t2- 
The reflection and transmission coefficients of this struc- 
ture must read 
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where r° and t° are the Fresnel coefficients of the bare 
ei|e2 interface, while 77 is the self-consistent amplitude of 
the wave scattered by the layer in response to both the ex- 
ternal field and the multiple reflections within the layer- 
interface cavity. The scattered field is obviously symmet- 
ric for a layer of thickness -C A, and it contributes to Eq. 
^ with direct (77) and interface-reflected (±r°?7) compo- 
nents, as well as the corresponding transmission (±t°r]). 
The absorbance A = 1 - |i?| 2 -Re{/}|T| 2 , obtained from 
the integral of the Poynting vector over planes far from 
the structure, is thus an analytical function of the com- 
plex variable 77, and its maximum admits the closed-form 
expressions 



and 
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for s andp polarizations, respectively, where / = (e 2 /ei — 
sin 2 8) 1 / 2 / cos 8 and 8 is the angle of incidence from 
medium 1. The dependence of the universal maximum 



absorption on 9 and 62/61 predicted by Eqs. ^ and 
Q is shown in Fig. [2]jb-c). Interestingly, it becomes 
100% under TIR from the high-index medium (i.e., for 
£i sin 2 8 > 62, leading to Re{/} = 0), indicating that it 
is possible to simply suppress the reflection when trans- 
mission is already prevented by TIR. 

Absorption in particle and disk arrays.- The univer- 
sal maximum is just an upper limit to the absorption, 
but its derivation does not provide any clue on how to 
structure the material to reach such a limit, or whether 
it is at all possible to obtain it. However, a systematic 
inspection of graphene nanodisk arrays reveals that the 
universal limit is indeed reached by engineering the disk 
size, doping, and spacing. We show several examples in 
Fig. [5Jd-e) , computed from a layer-KKR method [5D] , in 
which the disks are described through their multipolar 
polarizability obtained from BEM (solid curves). 

We can gain further insight into thin-layer absorption 
by describing each graphene nanodisk (or any other par- 
ticle) through the polarizability [Eq. Q]. For arrays of 
small period a compared to the wavelength, the reflection 
coefficient of the layer reduces to [2~Tj 
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where G is a lattice sum whose imaginary part is ex- 
actly given by Im{G} = S — 2(w/c) 3 /3, whereas its 
real part reduces to Re{G} « g/a 3 for A 3> a, with 
g = 5.52 (g = 4.52) for hexagonal (square) arrays. Here, 
S = 2iroj/cAcos9 (S = 2ttoj cos 8 / cA) for s-polarized (p- 
polarized) light and A is the unit-cell area. Equation |5| , 
together with Eq. ([T]), permits obtaining the transmission 
and reflection coefficients of the combined array-interface 
system analytically (they are simply given by Eq. ^ with 
77 = r(l ± r°)/(l — r°r)), leading to the dashed curves of 
Fig. |2jd-e), in excellent agreement with full simulations. 

In the symmetric configuration (ei = e 2 ), combining 
Eqs. ([T]) and jjjb, and noticing that k » c r , the maximum 
absorption (i.e., 50%, corresponding to r — Tl/2) is ob- 
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FIG. 2: Absorption by a layer of doped graphene nanodisks. (a) Scheme of the geometry under consideration. (b,c) Universal 
contours showing a maximum limit to the absorbance of a thin material layer sandwiched between two different media for s and 
p polarization as a function of incidence angle and dielectric contrast between both media. Light is incident from the ei medium. 
(d,e) Realization of maximum absorption with arrays of graphene disks [see (a)] of diameter D — 60 run, Ef = 0.4 eV, and 
different lattice spacings and angles of incidence (see labels). The lower-index medium has e = 1 in all cases. The numerical 
labels in (d,e) correspond to the positions of the symbols in (b,c). 



tained at a frequency u ~ uj p — 3gK r /4(u;a/c) 3 (slightly 
redshiftcd with respect to the single-disk resonance w p ), 
under the condition ( = 1/2, where we define the param- 
eter 



c 



.4 



cos 9, 



jCXt 
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s polarization, 
p polarization. 



This condition can be fulfilled by adjusting the spacing 
between disks, provided a^ax * s sufficiently large com- 
pared to the disk area. For asymmetric environments 
(f-i £2), the analysis becomes more involved, but simi- 
lar redshifts are obtained for large enough a^ax- 

The robustness of TLA under TIR is explored in Fig. 
|3j where the peak absorbance is represented as a function 
of incidence angle for various values of the lattice spacing 
a, using the same disks in all cases [see parameters in Fig. 
[3fb)]. Nearly TLA is observed within the 20°-45° range 
of incidence angles under critical coupling conditions [9] . 

Omnidirectional total light absorption. -The above pre- 
diction of TLA under TIR can be generalized to sys- 
tems in which the transmission channel is suppressed. In 
particular, we consider a graphene array placed above 
a metallic substrate. The graphene plasmons disappear 
close to a metal, but this problem is solved by introduc- 
ing an intermediate dielectric coating layer [see inset of 
Fig. |4|a)]. In this configuration, the metal can natu- 
rally act as a backgate element to dope the graphene. It 
is easy to see that reflection is canceled when the sub- 
strate (i.e., metal plus coating) reflection coefficient is 
r° = — r/(l ± 2r). The coating-layer thickness d can be 
adjusted to match the required phase in this expression, 
and the condition is simply stated as 1/|?~ | = |2±l/r| > 
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FIG. 3: Angular dependence of the peak absorbance for s 
(a) and p (b) polarization in arrays of graphene nanodisks 
sec Fig. g a)]. Each curve represents a fixed sample with 
parameters as shown by labels. 



1, where the equality applies to non-absorbing substrates 
(e.g., a dielectric under TIR, or noble metals in the 
IR). Figure ^b) shows TLA under this scheme for prop- 
erly adjusted values of a and d. Obviously, the phase- 
matching condition is periodically satisfied along the d 
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FIG. 4: Total absorption with graphene on a coated metallic substrate, (a) Normal-incidence absorbance by graphene disk 
arrays supported on a dielectric-coated gold surface (see inset) for various values of the dielectric film thickness (tdieiectric =2.1, 
£Au from |22|1. (b) Peak absorbance as a function of dielectric thickness d and lattice spacing a for the same disks as in (a). 
(c,d) Incidence-angle and photon-energy dependence of the absorbance under the conditions of (a) for d = 1.12 jitm. 



axis. Interestingly, the incidence-angle and polarization 
dependence of the absorption is very weak [Fig. |4jc-d)], 
pointing at the possibility of omnidirectional absorption. 
To understand this, we combine Eqs. ([I]) and ^ together 
with the relation |2 ± l/r| = 1 (i.e., the condition of no 
reflection for |r | = 1), from which we find two frequen- 
cies for TLA, lo — ujn — 3gK r /4:(uja/c) 3 ± {k/2)- s /1/Q - 1, 
under the condition C, < 1; only the right-most term de- 
pends on polarization and angle of incidence (through £) , 
so total absorption can be made omnidirectional when 
the second term in the resonant lo dominates (i.e., for 
(wa/c) 3 <C K r /n), provided c£f* x is of the order of A. 

Concluding remarks.- Planar textured materials can 
produce large light absorption, which reaches 100% when 
the transmission channel is suppressed. We predict that 
this effect takes place using a single sheet of doped pat- 
terned graphene either under TIR or when the carbon 
layer is deposited on dielectric-coated metal. Full absorp- 
tion in graphene opens interesting possibilities for light 
detection (e.g., via electron-hole separation [23], thermo- 
electric reading [23], etc.). This type of detector should 
inherit the electrical tunability of graphene, which can 
be used for direct, efficient spectral analysis of IR light 
with relative spectral resolution ~ n r / 1 n. Additionally, 
thermally heated samples can be useful as narrow-band 
IR sources in virtue of Kirchhoff 's law. We further pre- 
dict perfect absorption in planar particle arrays under 
the condition that the absorption cross-section of an in- 
dividual particle is comparable to the area of the unit 
cell. This condition is actually fulfilled by many non- 
ideal two-level atoms/molecules (e.g., NV centers). Our 
results could be tested with optically trapped atomic ar- 
rays, although graphene nanodisks are more convenient 
because they provide a practical implementation that is 
stable under ambient conditions. Noble-metal nanopar- 
ticle arrays can also exhibit similar effects in the near IR, 
with their spacing controlled through dielectric coating 



This work has been supported by the Spanish MICINN 
(MAT2010-14885 and Consolider NanoLight.es), Fun- 
dacio Cellex Barcelona, and the European Commis- 
sion (FP7-ICT-2009-4-248909-LIMA and FP7-ICT-2009- 
4-248855-N4E). 



[1] 

[2 
[3 

[4: 

[6 
[7 



[9 
[10 

[H 
[12 
[13 
[14 
[15 

[16 
[17 

[18 

[19 
[20 

[21 

[22 



6044 



Corresponding author: J.G.deAbajo@cs ic.es| 
K. Mizuno et al, Proc. Natl. Academ. Sci. 106, 
(2009). 

P. B. Clapham and M. C. Hutley, Nature 24, 281 (1973). 
M. C. Hutley and D. Maystre, Opt. Commun. 19, 431 
(1976). 

M. Neviere et al, in Proc. of ICO- 11 Conference, Madrid, 

Spam (1978), pp. 609-612. 

E. Popov et al, Opt. Express 16, 6146 (2008). 

E. Popov, L. Tsonev, and D. Maystre, Appl. Opt. 33, 
5214 (1994). 

Y. P. Bliokh, J. Felsteiner, and Y. Z. Slutsker, Phys. Rev. 
Lett. 95, 165003 (2005). 

N. I. Landy et al, Phys. Rev. Lett. 100, 207402 (2008). 
T. V. Teperik et al, Nat. Photonics 2, 299 (2008). 
K. S. Novoselov et al, Science 306, 666 (2004). 
K. S. Novoselov et al, Nature 438, 197 (2005). 

F. Bonaccorso et al, Nat. Photon. 4, 611 (2010). 

C. F. Chen et al, Nature 471, 617 (2011). 

K. F. Mak et al, Phys. Rev. Lett. 101, 196405 (2008). 
M. Jablan, H. Buljan, and M. Soljacic, Phys. Rev. B 80, 
245435 (2009). 

F. H. L. Koppens, D. E. Chang, and F. J. Garcia de 
Abajo, |arXiv:1104.'2068V L 

F. J. Garcia de Abajo and A. Howie, Phys. Rev. B 65, 
115418 (2002). 

D. V. van Coevorden et al, Phys. Rev. Lett. 77, 2412 
(1996). 

B. Wunsch et al, New J. Phys. 8, 318 (2006). 
N. Stefanou, V. Yannopapas, and A. Modinos, Corn- 
put. Phys. Commun. 113, 49 (1998); 132, 189 (2000). 
F. J. Garcia de Abajo, Rev. Mod. Phys. 79, 1267 (2007). 
P. B. Johnson and R. W. Christy, Phys. Rev. B 6, 4370 
(1972). 



5 



[23] M. C. Lemme et al, arXiv: 1 012.4745vl. Chem. B 105, 3441 (2001). 

[24] J. C. W. Song et al., |arXiv:1105. 1142^ 1. 

[25] T. Ung, L. M. Liz-Marzan, and P. Mulvaney, J. Phys. 



